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We investigate the modification in mesoscopic electronic transport due to electron-electron inter- 
actions making use of scattering states. We demonstrate that for a specific (finite range) interaction 
kernel, the knowledge of the scattering matrix is sufficient to take interaction effects into account. 
We calculate perturbatively the corrections to the current and current-current correlator; in agree- 
ment with previous work, we find that, in linear response, interaction effects can be accounted for by 
an effective (renormalized) transmission probability. Beyond linear response, simple renormalization 
of scattering coefficients is not sufficient to describe the current-current correlator, as additional cor- 
rections arise due to irreducible two-particle processes. Furthermore, we find that the correlations 
between opposite-spin currents induced by interaction are enhanced for an asymmetric scatterer, 
generating a nonzero result already to lowest order in the interaction. 
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I. INTRODUCTION 

Transport through mesoscopic systems in the presence 
of particle scattering and interaction has been at the fo- 
cus of experimental and theoretical studies during sev- 
eral decades. The scattering matrix theory reducing the 
problem to asymptotic states^^ proved to be an efficient 
tool characterizing transport in the presence of single- 
particle scattering and for non-interacting systems. Be- 
yond calculation of t he current, the scattering matrix for- 
malism was usecP^ to determine the noise properties of 
devices and provided the complete probability distribu- 
tion of transmitted charge through the calculation of the 
full counting statistics (FCSpEI. Several reviews dealing 
with the scattering matrix approach to electron transport 
are available today^H 

Effects of electron-electron interactions give rise 
to numerous interesting phenomena, e.g., Coulomb 
blockade^, Kondo effectPSI, or the 0.7 anomaly^ in trans- 
port, or entanglement and decoherence within the con- 
text of quantum processing in mesoscopic systems^!. In- 
teractions arise locally and thus involve the particle wave 
function within the domain where interactions take place. 
A natural way to tackle the combined problem of scat- 
tering and interactions is to start with the interacting 
problem (e.g., the quantum dot) and include scattering 
thereafter, as in the tunneling approach which treats the 
scattering part perturbatively. In the present paper, we 
take an alternative approach which describes the scat- 
tering part exactly for any scatterer. We then have to 
compromise with the interaction, for which we choose a 
particular kernel which is constant over the scattering 
region. This scheme allows us to express the full inter- 
acting Hamiltonian through the exact scattering coeffi- 
cients of the non-interacting system. The approach can 
be extended to the case of electron-photon/plasmon and 
electron-phonon interactions, where the wave function 
of the scattered low-energy mode (rather than the two- 
particle interaction kernel) is approximately constant 



over the region of scattering and interaction. This idea 
has been exploited before in the study of photon/plasmon 
generation by electron scattering, as discussed in Refs. 

Accounting for interaction effects in electronic trans- 
port is a highly non-trivial task and general results are 
sparse, with notable exceptions such as the generalized 
Landauer formula due to Meir and WingreerPS however. 
Still, considerable progress can be made for specific sys- 
tems: Studies of the (AndersorP^) impurity model for 
quantum dots provide non-perturbative results for the 
non-linear I-V characteristics^, for the shot noise^, and 
perturbative results (and beyond) are available for the 
full counting statisti cs^. The same quantities, current- 
voltage characteristic! 2 ^!], shot-noise^, and full counting 
statistics^ , have been found for the interacting resonant 
level mode l 25 * 26 ! The more general case of a quantum 
point contact with an arbitrary two-particle interaction 
has been analyzed with the help of field-theoretic meth- 
ods: the current and noise have been determined by Gol- 
ubev et alW^ within a saddle-point approximation and 
the renormalization group flow of the transmission coef- 
ficients has been found by Kindermann and Nazarov^ 
in their study of weak interaction effects on full count- 
ing statistics. In our approach, we express the interac- 
tion through the non-interacting scattering coefficients 
of an arbitrary mesoscopic conductor, yet restricted to a 
capacitive type of Coulomb interaction. Providing this 
starting point, this interacting Hamiltonian can be used 
in a simple perturbative calculation as done below, but 
could also be fed to a more sophisticated technique as 
those developed in Refs. 1271281 

Much interest has been devoted to the question 
whether the Landauer formula for the current survives in 
the presence of interactions, a question that our approach 
is quite suitable to address. For systems where the inter- 
action among electrons is restricted to a finite region of 
space, Meir and WingreerP^ developed a generalization of 
the Landauer formula for the current, which was recently 
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shown to reduce to the form of the Landauer formula 
via introduction of an effective transmission^. Further- 
more, the analysis has been extended to include effects 
of interactions in the leads®'. Gogolin and KomnilPSl, 
devising a method providing the generating function of 
the full counting statistics, have found that, within lin- 
ear response, effects of interaction in the Anderson model 
can be accounted for with a renormalization of scatter- 
ing coefficients. However, it is still under debate whether 
all effects of interaction can be included via an effective 
transmission probability in the Landauer formula®! 

Below, we consider two (semi-infinite) non-interacting 
leads connected via a finite region of space where single- 
electron scattering and interaction between electrons 
takes place. We focus on an interaction between elec- 
trons which is restricted to a finite region in space and 
furthermore constant in the region of single-electron scat- 
tering. In this case, the matrix elements in the interac- 
tion Hamiltonian can be calculated from the knowledge 
of the asymptotic behavior of the single-particle wave- 
functions (scattering matrix) . Our study is based on the 
Keldysh formalism^ anc j includes the non-equilibrium 
situation with a finite applied bias voltage V. We treat 
the scattering problem exactly by writing the initial den- 
sity matrix p(to) (with to — > — oo) in terms of the scat- 
tering states of the non-interacting system. Assuming an 
arbitrary scatterer, we limit ourselves to a perturbative 
expansion in the interaction and stay with the lowest or- 
der (in which case the introduction of Green's function as 
done below is not a must). Adopting a specific scatterer, 
e.g., a single resonance level, we push the perturbative 
expansion to higher order with a selected choice of dia- 
grams. 

In a first application of our approach, we determine the 
interaction-induced correction to transport and noise: to 
lowest order in the interaction strength Uq and to all or- 
ders in the applied bias V, the corrections to the single- 
particle Green's function G and the transport current 
j can be described through renormalized scattering am- 
plitudes. On the other hand, the corrections to the two- 
particle Green's function and the current-noise cannot be 
cast into the form of renormalized scattering coefficients 
in general. This is due to the presence of irreducible 
two-particle processes (vertex corrections). Calculating 
the current and noise from the single- and two-particle 
Green's function, we find that in the linear voltage regime 
V all the effect of interaction can be accounted for by an 
effective transmission (to order Uq), in agreement with 
previous findings by Gogolin and KomnilP2. In the non- 
linear voltage regime, this is not the case any more and 
interaction effects beyond a simple renormalization of the 
scattering properties appear. Furthermore, it turns out, 
that the corrections to the two-particle Green's function 
generate a correlation between opposite-spin currents in 
an asymmetric device which appears to order UqV 2 7 while 
in a symmetric system®! sucn a correlation arises only in 
order UqV 2 . Again, our approach turns out to be quite 
suitable in unravelling such additional correlations due to 



an asymmetry of the scatterer, a question that has not 
been addressed before. 

The paper is organized as follows: In Sec. [XXJ we 
introduce the Hamiltonian, briefly explain the non- 
equilibrium formalism used in the paper, and demon- 
strate how the interaction matrix elements simplify due 
to our particular choice of the interaction kernel. In 



Sec. Ill we calculate the lowest order (in interaction) 



correction to the current and noise correlator and find 
the interaction-induced correlation between opposite- 
spin currents. We then evaluate the results for the spe- 
cific case of a single resonance level in Sec. |IV[ including 
a resummation of diagrams providing the mean-field cor- 
rection to the single-particle Green's function. A sum- 
mary and conclusions are given in Sec. |Vj 



II. FORMALISM 

Our one-dimensional system involves two half-infinite 
leads connected by a central region where particles are 
scattered and is described by the Hamiltonian 



H — H ] 



kin 



v, 



(1) 



where the potential V(x) gives rise to single-particle scat- 
tering. Making use of the Lippmann-Schwinger (LS) 
scattering states | we treat the scattering prob- 

lem exactly; the index a = L/R distinguishes left/right- 
incoming states and we choose the wavevector k > 0. 
The LS scattering states are eigenstates of the Hamilto- 
nian H n , 



(H kin + V)\(p ak ) = E k \cp ak ), 



(2) 



with energy E k . In the following, we will consider 
a spectrum linearized around E F = hu F , i.e., E k = 
E F + hv F (k—k F ) (this is a convenient simplification rather 
than a necessity). The asymptotics of the LS states 
(ip a k(x) = (x\ip ak )), 



e ikx _|_ gifcdg ikx^ x _^ 

t k e lkd e lkx , x^ 

t k e lkd e~ lkx , x — > 

e~ lkx + rR k e lkd e tkx , x —> 



-oo, 

00, 

-oo, 

00, 



(3) 



is determined by the scattering matrix 



Sh — 



( Thk t k 

V tk r Rk 



(4) 



where tjj. — tji k = t k in our time-reversal symmetric 
system. Different from more standard formulations, we 
have defined the scattering coefficients with additional 
phase factors exp(ifcc?) as this will simplify our expres- 
sions later. The Lippmann-Schwinger states satisfy the 
orthogonality 



{<Pak\¥>j3q} = 27T(5(fc - q)5 a p 



(5) 
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and, in the absence of bound states as assumed here, the 
condition of completeness. The Hamiltonian H Q then can 
be rewritten in terms of the LS states as 



Efdk * 
I — Eu c , c , 
/ 2yj- ctko aka 



(6) 



with the creation (annihilation) operator cj, (c , ) of 
LS scattering states with spin cr. With a voltage U ap- 
plied to the system, we assume a non-equilibrium steady 
state at an initial time t (t — > — oo) which we describe 
by the initial density matrix 



Pa = p(t ) = exp 



£ 

Oi.CT 



dk 
2^ 



E k 



(7) 



with the chemical potential (J-l/r — E F ± eV/2 in the 
left/right lead. 



A. Interaction matrix elements 

We consider an interaction between electrons described 
by the Hamiltonian 



2 ^ 



(8) 



with U(x,x') the interaction kernel. Here, we focus on a 
situation where the interaction can be assumed constant 
within the region of single-particle scattering. Typical 
physical systems we have in mind are quantum dots (with 
constant interaction inside the dot and vanishing outside) 
and quantum point contacts (where the interaction can 
be assumed large and constant between electrons in the 
quantum point contact region). Inspired by the capaci- 
tive form Q 2 /2C (with the charge Q and capacitance C) 
we choose an electron-electron interaction kernel 



U(x,x') = U f(x)f{x'). 



(9) 



Using the Lippmann-Schwinger scattering states, we ex- 
press the field operators as 



$<r{x) = j 7j- (f a k(x)c aka 



and can rewrite the Hamiltonian as 

£j U ^-^ ^-^ f dk dk 1 dq dq' 



(10) 



2 ^ ^ I 2tt 2tt 2tt 2tt 

era' aa' ' (5(3' ' 
x ^-a'k' ,ak^[3'q' .,(3q ^a_'k' a^' q' a'^^qa'^Oika 



(ii) 



with the matrix elements 



-4, 



a k,p q = / dxf(x)Lp* ak (x)y pq {x). 



(12) 



Assuming a scattering potential V(x) which vanishes out- 
side the region [— d/2, d/2], we consider a specific kernel 
with 

f(x) = e(d/2-\x\), (13) 
allowing us to rewrite the matrix elements in Eq. (fl2|) as 



A a k,/3q = (<Pak\<PPq) 



dx V* a h{x)<pp q {x) (14) 



with C d = K\ [-d/2, d/2}. The first term is given by the 
orthogonality condition Eq. ^ of the LS states, while 
the second term can be calculated using the asymptotic 
form Eq. ^ of the LS states. As a result, the overlap 
matrix elements can be expressed through the scattering 
coefficients alone, 



a _ ■ 1 ( r ak r aq + ^Aj) j(k-q)d/2 

-<*-ak.aq — ^ , & ; 

k — q 



k — q 



(15) 



where we have dropped terms suppressed by the small 
parameter l/(dk F ) and a ^ 13. For coinciding wave vec- 
tors k = q the result reads (again, a ^ /3) 



A a k,ak — ~i{ r *ak^k r ak + Cfc^fc*afc)> 

A ak ,pk = -i{r* ak d k tg k +t* ak d k ra k ). 



(16) 



It is straightforward to generalize this scheme to situa- 
tions where the kernel f(x E [—d/2, d/2}) — 1 decays 
smoothly outside the (scattering) re gion [—d/2, d/2]. For 
the interacting resonant level modeP^2£l, the interaction 
can be described by the Hamiltonian Eq. ^ with a ker- 
nel U(x,x') — Uog(x)f(x') and the function g(x) finite 
in the regions nearby but outside the dot; the interac- 
tion Hamiltonian Eq. (11) then can be expressed by the 
product of matrix elements A ak ^ q and a similar quantity 



with f(x) replaced by g(x) in Eq. (12). In the following, 
we concentrate on the capacitive interaction Eq. ^ with 
the specific choice Eq. ( [l3| . 

Another situation where the knowledge of the scatter- 
ing coefficients suffices to evaluate matrix elements is the 
interaction with long wave-length bosonic modes. As a 
specific example, we consider plasmonic excitations of a 
two-dimensional electron gas (2DEG) induced by scat- 
tering electrons, e.g., at a point contact; the interaction 
is described by the HamiltoniarP^, 



-ihd x 



A(x) Mx), (17) 



where A(x) is the vector potential of the electromagnetic 
field induced by the two-dimensional plasmons. For the 
lowest transverse component of the plasmon wave func- 
tion <pk(x) one can write, 



A(x) 



£ 

k 



k^y 



/27rfic 2 \i/2 



'{bMx)-b\4>l{x)}, (18) 
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where b k (b k ) are bosonic creation (annihilation) opera- 
tors for the plasmon modes with wave vector k and linear 
dispersion uj k = v p \k, with the plasmon velocity v p \ and 
a dimensionless geometrical factor 7 ~ 1. 

Typically, v p \ 3> v F , and the conservation of energy 
5E = Vpihkpi 3> v F hk implies that the plasmon wave 
length is much larger than the Fermi wave length, A p i S> 
A F . The size of the plasmon then exceeds, by several or- 
ders of magnitude, the scattering region (the size of the 
quantum point contact) and thus the interaction between 
the plasmons and electrons is dominated by the asymp- 
totic scattering region where the electronic wave function 
has the form ([3|. Neglecting the spatial dependence of 
the plasmon wave function 4>(x) = const., one again can 
express the interaction Hamiltonian (17 1 solely in terms 



of the non-interacting scattering matrix. 



B. Correction to Green's functions and scattering 
matrix 

We are interested in calculating the transport prop- 
erties of a device, such as the mean current and the 
current-noise correlator. Within the Green's function ap- 
proach, we need to know its single-particle version for 
the calculation of the current, while for the calculation 
of the current-current correlator, knowledge of the two- 
particle Green's function is required. Below, we first de- 
rive the interaction correction to the single-particle and 
two-particle Green's functions, before focusing on the 
transport properties in the next section. 



forrrPl 



W^afc>')M- (20) 

dk <p ak {x)tp* ak (x') 



with the wavevector k u 



2ixi k — fc w — iS 
k F + (id — lu f )/v f , the matrix 



-n a (io) -n a (oj) 
1 — n a (ut) 1 — n a (oj) I ' 



(21) 



and the Pauli matrix a z . Here, n L/a (uj) is the Fermi 
distribution for the left/right lead. The integral in the 
expression Eq. ( p0| of the non-interacting Green's func- 
tions can be explicitly calculated for coordinates outside 
of the scattering region, i.e., |x|,|x'| > d/2, by insert- 
ing the asymptotic wave functions Eq. ([3j). E.g., for 
x > x' > d/2, we obtain 



G^(x,x',u>) = -^e lk ^ x - x '\M L ^+a z ) 



(22) 



_| l e -ik w (x-x') £j 



V F 

—i 



Ru 



'-R k y k ^ x ~ x '\M J 



Ru 



Mr 



' Rk, 



' Rk, 



a ik^{x+x' — d) I 



1. Single-particle Green's function 

The single-particle Keldysh Green's functions^! are de- 
fined as 



.(j (x , t, x . t ) 

= -i(T K [$v(x,U)'ipl(x' ,t' u ,)S K _ 



(19) 



where the expectation value of an operator A is defined 
as (A) = Ti{ Ap } with the trace taken over all many- 
particle states of the system and p is given by Eq. Q. 
We will skip the spin-index for the single-particle Green's 
functions in the following, as these are equal for both 
spin directions. The time-dependence of the field oper- 
ator ip(x,t) in the interaction representation is governed 
by the free Hamiltonian Eq. u (z/) is a Keldysh in- 
dex ± telling on which branch of the Keldysh contour 
C = (t a ,oo) + U (00, t )— the time argument t (t') is lo- 
cated (with t —00). Tk is the Keldysh time-ordering 
operator and Sk is the Keldysh evolution operator in the 
interaction representation. 

The non-interacting Green's functions in frequency 
representation G°l, (x,x',uj) can be written in matrix 



I 2 and T k ^ — \t Lku _\ 2 and we have 



where R k ^ = \r RK 
used T k ^ + R k = 1 to simplify the expression. Eq. ( 22 ) 
expresses the non-interacting Green's functions outside 
the scattering region in terms of scattering matrix coeffi- 
cients; as we will see below, the (lowest order) corrections 
due to interaction assume the same structure and can be 
formulated in terms of a renormalization of the scattering 
coefficients. 



v 
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FIG. 1: The Hartree and Fock diagrams give rise to the first 
order corrections to the Green's functions given by Eq. ( 23 1 
and plk. 



The first-order correction G m to the single-particle 
Green's function is given by the first-order expansion of 
Eq. ( |19[ ) in the interaction Hamiltonian. There are two 
contributions in lowest-order, a Hartree and a Fock term, 
cf. Fig. [T] for which we find the expressions 
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G^(x,x',ij) = -j j dydy'U(y,y') f ^ &°>(x, y, u)* z G™ (y', y',Cj)G™(y, x',u), (23) 



G^(x, x',cj)= t - f dydy>U(y,yi) / g &°>(x, y, u)& s G^_{y, y' ,u>)&°\y' , x',u). 



(24) 



r 



The real-space non-interacting Green's functions 
G m (x, x' , lu) can be expressed through the scattering 
states using 

G<°>0r,s»=E / ^SSMVafctoVafcOO, (25) 



simpler results 



(ih) _ U N{ n t 

Tak ~ h^ dkrak ' 



(28) 



where are the non-interacting Green's functions 

(in fc-space) 



ffikM = HkboXMw + a z ) - g£{u)M Q 



(26) 



ft I A 

with the retarded/advanced Green's functions g k Joj) = 
(uj— ui k ±iS)~ 1 and uj k = E k /h. For our kernel Eq. (B| and 
the specific choice of f(x) in Eq. (13), the integrations 
over the internal coordinates reduce to integrations of two 
scattering states within the region [— d/2,d/2]. We re- 
place these integrals by integrations over the asymptotic 
regions as in Sec. |II A] before in order to avoid the integra- 
tion over the scattering region. As a result, these overlaps 
can be expressed through the matrix elements A ak ^ q de- 
fined in Eq. (15). The expressions for the Hartree and 



Fock corrections can be evaluated for all coordinates out- 
side of the scattering region. Comparing for each choice 
of coordinate positions (to the left /right of the scatterer) 
the corrections with the non-interacting Green's function 
itself, we recognize that we can cast all effects of interac- 
tion into a renormalization of the scattering coefficients. 
E.g., for x,x' > d/2, the corrections assume the form 
of the three last terms in Eq. (22) and we can define 



the renormalization r Rk to the reflection coefficient r Rk 



in a consistent manner such as to capture all three cor- 
rections. The corrections to the scattering matrix again 



involve Hartree- and Fock contributions, t^l 



4(1H) , |(1F) 

l ak 



and r (1 ? 

ak 



' ak 



' ak ' 



the Hartree terms are given by 



Lk ~ 



2Uo 



Rk 



2Uo 



—npqiitk A Lk , Lk Ap qtl3q (27) 
P 71 

+ irRk A Rk ,LkAp q ,f3 q ) , 

^ / Y"np q (ir Rk A Rk , Rk Ap q fi q 
p J n 

+ itf- ALk t MkAp qt p q ) , 



and the corresponding corrections tp k and 



Lk 



are ob- 



tained by interchanging R o L. The q- and fc-dependent 
parts in the above expressions factorize and we find the 



with the particle number iVj nt in the interaction region 
given by 

d/2 

N int = J2 J dy{i>t(y,t)i> a (y,t)) 

" -d/2 

— (n Lq A Lq x q + n Rq A Rq . Rq ) , (29) 



where n ak = n a (E k ) and the factor 2 accounts for the 
two spin directions (note that the fc-dependent factors in 
Eq. (27) reduce to a derivative of the scattering ampli- 
tudes, i.e., the brackets in Eq. (27) reduce to Ap q ^ q d k tk 
and Af3 qt p q d k r Rk , respectively). 

Given the corrections to the scattering amplitudes, we 
can calculate the corrections to the transmission and re- 
flection probabilities T^ 1H) 



and R[ 1H) 



fp(lH) 
~^k ' 



(30) 



The result is easily interpreted: the Hartree interaction 
effectively shifts the scattering potential with respect to 
the energy of the incoming particle and thus the trans- 
mission characteristics of the scatterer is shifted accord- 



ingly; Eq. ( 30 1 then describes the lowest order correction 
of T k due to a shift in wavelength by Sk = —U N[ nt /Hv F . 

The Fock correction to the scattering matrix assumes 
the form 



Lk 



Rk 



Hi 



J ^np q (itk ALk,p q Ap q>L k 



(31) 



ir Rk A R k,pqAp qi Lk) , 



^ / -^- n Pq{ ir RkA Rk ,p q Ap q . 



2tt 



Rk 



it k A 



Lk,f3qAp qt Rk) 



The corrections i^j? and can be found by exchanging 
R <-> L in the expressions above. While for the Hartree 
term the corrections to the transmission coefficients are 
equal, i.e., t 1 ^ = t^f , this is in general not the case for 



G 



the Fock contribution (however, equality holds in equilib- 
rium) . The corrections to the transmission and reflection 
probabilities T^ 1F) = -R ( k 1F) are 



2tt (k - q) 2 



(32) 



The modifications to the current and the current-current 
correlator due to the renormalization of the single- 
particle Green's functions will be discussed in section |III| 
below. Higher order perturbation theory might gener- 
ate terms of different form which cannot be incorporated 
within a simple renormalization of the scattering coeffi- 
cients. 



with the single-particle Green's function G(l, 1') 



G^^ixi^ti, x[,t[), cf. Eq. (19 1, and similar for the 



other coordinates (again, we ignore the spin index). The 
non-interacting reducible two-particle Green's functions 
{Ga^ d) ) {0) an d corrections due to interaction (e.g., 
(Q^ d) ) (1} ) can be calculated from the non-interacting 
single-particle Green's functions G (0) and the correction 
G (1) to the single-particle Green's function. 





as 



2. Two-particle Green's function 
The two-particle Keldysh Green's functions are defined 

CT7r (l,2;l',2') (33) 
= (-*) a (75 f (^(l)^(2)^(2')^(l')^ic)) ) 



where the variables i — 1,2,1', 2' represent both space 
and time coordinates as well as the Keldysh index, i.e., 
i = (xi,ti,i>i). The two-particle Green's functions can 
be split into a reducible and an irreducible part, Q an = 
GaZ d) + Go* 5 > where the reducible part can be expressed 
through single-particle Green's functions as 

er , (l,2;l',2') = G(l,l')G(2,2') (34) 
-<UG(1,2')G(2,1') 



FIG. 2: There are two processes contributing to the correc- 
tion to the irreducible two-particle Green's function in first 
order given by Eq. ( 35 1. While the second process only con- 



tributes for equal spins, the first one also gives a contribution 
for opposite spins. 

The remaining irreducible part cannot be ex- 

pressed through single-particle Green's functions and is 
given by the connected diagrams of Eq. ( [33] ) . The lowest- 
order contribution to the irreducible two-particle Green's 
function involves two contributions, cf. Fig.[2j and can be 
evaluated for coordinates outside the scattering region 
and arbitrary Keldysh indices using the non-interacting 
Green's function in Eq. |20[ ); to first order in the inter- 
action (note that = 0)> we obtain in frequency 
representation with i — [x^uii,Vi), 



(^) (1, (1,2;1',2') 



.U 



d/2 



i^r2ir6(uJi + ui2 — u[ — ui' 2 ) J dy dy' 

-d/2 



E 7 [g ( ^(xi, y, ^)G^(x 2 ,y', ^ 2 )G<°4 (y>, ^^)G™, (y, ^M) 

7 

- S^G^xuy, wi)G<« (za, y', w 2 )G™ ; (y, x' 2 , 4)G^ ; (y', ^M) 



(35) 



The integrals over the region [— d/2, d/2] are calculated 
in the same way as in Sec. |II A| above. 



III. CORRECTION TO CURRENT AND NOISE 

Next, we determine the interaction corrections to the 
current and current-current correlator as calculated from 
single-particle and two-particle Green's functions; the ir- 
reducible contributions to the latter will give rise to spin 
correlations in an asymmetric device. 
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A. Current 

The current j(x,t) derives from the single-particle 
Green's function via 



j{x,t) = 2t> x G + 4x,t,x',t), 



(36) 



where T> x = (eh/2m)(d x > — d x ) and the limit x' — > x 
is taken in the end. The factor 2 accounts for the two 
spin orientations of the electrons. We determine the cur- 
rent to the right of the scattering and interaction region, 
x > d/2. As the effect of interactions on the single- 
particle Green's function can be expressed through a 
renormalization of the scattering coefficients, the same 
applies to the current and we arrive at the well-known 
Landauer formula 

jV) +j W=2evJ ^ [Tk +T^}(n Lk -n Rk ) (37) 



with T k = T£ M> +T£*> given in Eqs. (30 1 and (32). 



B. Current-current correlator 

Let us now turn our attention to the current-current 
correlator 

S , - 7r (w,xi,x 2 ) (38) 
= J dr((T K (Uxi,T-)Ux2,0 + )S K )))e iUT 

behind the scatterer, x\,x 2 > d/2. Below, we focus 
on the zero-frequency correlator which is independent 
of X\ and X2, i.e., S a7T (0) = S aw (0, X\, X2), as explained 
in Ref. 1331 The current-current correlator can be ex- 
pressed in terms of two-particle Green's functions and 
thus can be split into reducible and irreducible parts, 
S aiT = S ( J° d) + S'^ } . The reducible part vanishes for op- 
posite spins a 7^ tt, while for equal spins it is expressible 
in terms of single-particle Green's functions, 

dui 
— ( 
2tt 

(39) 

The irreducible part contributies to both, equal and op- 
posite spins, and can be expressed as 



S^(0) = -V XX V X2 I ^-G_ + (x 1i x' 2 ,lo)G + _(x 2 ,x' 1 ,u). 



S™(0)=T> xl V t 



dujx duJ2 (itr) 
(2tt) 2 CT7r v " ' 



2 + ,l'_,2' + ), (40) 



where we explicitely indicated the Keldysh indices and 

LJx> t 2' = ^1,2- 

In the non-interacting case, the irreducible Green's 
function vanishes and the current -cu rrent correlator is 
given by the well-known expression^ 

S^(0)=5^e 2 v F [ g (41) 



T k (n Lk (l - n Lk ) + n Rk (l - n Rk )) 
+ R k T k {n Lk {\ - n Rk ) + n Rk (l - n Lk )) 



There are two kinds of (first-order) interaction correc- 
tions to the noise correlator: those due to the reducible 
part S^ d) can be expressed through the renormalization 
of the transmission probability T k in Eq. (41 ). The con- 



tribution to the irreducible part gives rise to addi- 
tional correlations which go beyond a simple renormal- 
ization. 



Reducible corrections 



The expansion of the reducible part Eq. (39) in the 



interaction produces corrections originating from those 
in the single-particle Green's functions and we find 



dk 
2^ 



(42) 



2T k T { k ] (n Lk (l - n Lk ) + n Rk (l - n Rk )) 
+ (1 - 2T k )T { k 1] (n Lk (l - n Rk ) + n Rk (l - n Lk )) 



which corresponds to the non-interacting expression 
Eq. (41) with a rcnormalizcd transmission probability 
T k +W } 



: k and expanded to first order in T k 



2. Irreducible corrections 

Restricting ourselves to zero-temperature ($ = 0), we 
find the first-order contribution to the irreducible two- 



particle Green's function Eq. ( 35 ) with all coordinates 



to the right of the barrier x j , x\ > d/2 and evaluate the 
expression Eq. ( 40 ) for the irreducible noise correlator to 
arrive at 



C(irr)(l) 
(T7T 



(°)Lo 



2^0 . / N 

-e -r-sign(/i L - /Ur) 



dki dk 2 



irlk! *kj r Lk 2 tk 2 A Rk-i ,Lki A R k 2 Xk 2 
- ^Rkih^Rk^h A Lk u Rk^ A Lk 2 ,Rk^j 

- 5a7r(ir* Lkl t kl r* Lk2 t k2 A RkltLk2 A Rk . 2 ^ Lkl 
- ^Rfe 1 *fe 1 r Hfe 2 *fe 2 A LklyRk2 A Lk2tRkl 



(43) 



where A: max/min = max/min(fc F]j /A; pa ). The above correc- 
tion vanishes for the symmetric setup where r^ k — r Rk 
and A Rkt L q — AL k ^ Rq ; on the other hand, interactions 
combined with an asymmetric scatterer, a case that has 
not been addressed previously, generate a finite correla- 
tion already in lowest order in Uq. 

In summary, for equal spins, the correlation (or noise) , 
due to scattering and fermionic statistics already present 
without interaction, is modified due to interactions via 
the renormalization of the transmission (in the reducible 
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FIG. 3: (a) A (spin-up) electron on the dot shifts the reso- 
nance level for incoming (spin-down) electrons, thereby mod- 
ulating the (spin-down) current across the scatterer. (b) For 
an asymmetric scatterer, the escape probability of the (spin- 
up) electron to the left (Pl) and right (Pr) leads are different, 
generating a net (spin-up) current which is correlated with the 
(spin-down) current modulation. 



contribution) and by the additional irreducible contribu- 
tion given by Eq. ( [43] ) (note that this first-order contribu- 
tion may vanish for equal spins due to the Pauli principle, 
as, e.g., for a single-level quantum dot). Opposite-spin 
correlations are absent in the non-interacting case, and a 
finite result is generated to first order in the interaction 
for an asymmetric scatterer by the first pair of terms in 



Eq. (43) and to second order in Uo for a symmetric scat- 
terer. We attribute the origin of the asymmetry-induced 
correlations to the following effect (we restrict the dis- 
cussion to the case of an asymmetric quantum dot, cf. 
Fig. [3|: First of all, we note that, by virtue of Pauli's 
exclusion principle, all interaction-induced modifications 
of the noise correlator are due to the presence of spin 
a particles on the dot modifying the current of spin —a 
particles traversing the dot. The finite correlation arises 
from the change in current (5J_ CT due to the presence of 
particles 5Ni ntlJ on the dot, cf. Fig. [3] (a). For a symmet- 
ric scatterer, the decay of the charge 5Ni ntr7 to the un- 
blocked scattering state is symmetric and the undirected 
current leads to a vanishing current-current correlator. 
On the other hand, in an asymmetric dot, the charge 
SNinta decays in an asymmetric manner and its current 
8I a is correlated with 8I- a , cf. Fig. |3](b). 



IV. SINGLE RESONANCE LEVEL 

The above results, which apply to an arbitrary scat- 
terer, can be simplified and further developed when 
choosing a particular scatterer. Here, we consider a sim- 
ple example in the form of a scatterer with a single reso- 
nance level at energy E TCS = E F + Hv F (k TOS — fc F ), located 
at x — with vanishing extent d — > 0. This single-level 
quantum dot is described by the scattering matrix Sk 
with amplitudes 



tk — 



rLk 



TRk 



Ak + 77 
Ak + i?/7 
Ak + «7 
Ak — ir/j 



(44) 



Afc + 17 ' 

where the wavevector Afc = k — fc rca is measured rel- 
ative to the resonance level; the level width hv F j is 



(b) 



E«( W )= Q 



+ 



FIG. 4: (a) The two-particle vertex T^/ 2/ „, is anti- 
symmetrized with respect to the two incoming legs and to 
the two outgoing legs, (b) The first-order self-energy T,^(cu) 
consists of a Hartree and a Fock contribution. 



parametrized by 7 and the dimcnsionless parameter 77 G 
[—1,1] characterizes the asymmetry of the system. This 
scattering matrix generates a transmission probability 

T fc H**l a = (W)/(l + (A*/7) 2 )- 

The scattering matrix of the single resonance level 

model gives rise to matrix elements A a k t p q of a particu- 
larly simple form, i.e., 



A ak ,p q = a* a 4>* k ap(f> q , 



(45) 



with cll/r — =j= V and 



interaction Hamiltonian Eq. (11) then factorizes in the 



% /7/(Afc + i 7 ). The 



momenta, with individual terms given by 



(46) 



x4 



ai<j' 1 °<j 2 a! 2 C a' 1 k' 1 <T' 1 C a' 2 k' 2 a- 2 C a 2 k 2 a 2 C a 1 k 1 a 1 I 

and integration, summation over all wave vectors, lead, 
and spin indices providing the full interaction Hamilto- 
nian H m t • Introducing the anti-symmetrized two-particle 
vertex 

r i'2',2i = ^oa^a^a^a,,^^^^^^^,^^ (47) 

with i = 1, 2, 1', 2' representing the multi-index 

(oii,ki,ai) and 



S. 



(T 1 a 2 )0"2<7 1 



(48) 



the individual terms in Eq. ( 46 1 assume the form 
1 



'lint — ^l'l' : 21^ a ' 1 k' 1 a' 1 ^a 2 k 2 a 2 ^a 2 k 2 a 2 ^a 1 k 1 a 1 - (49) 

The Hartree and Fock digrams in Fig. [T] then collapse 
into one diagram, cf. Fig. [4j 

Given the simple structure of the vertex ( [47) ) , we can 
easily sum the Hartree-Fock diagrams for the single- 
particle Green's function (cf. Fig. [4^b)) with the help 
of the Dyson equation 



at 



^ J 2tt 2tt 



(50) 



9 



The first-order (Hartree-Fock) contribution to the self- 
energy X^ 1 ) assumes the following form (we use the 
Keldysh matrix notation and suppress the spin-indices 
as the self-energy is oc S a ' a ) 



(51) 



with U = UoNi n t/2 and N int , cf. Eq. (29), is the number 
of electrons on the level, 



iV int = 



d{q/l) (1 ~ v) n Lg + (1 + ?7)n flg 



2tt 



1 + (Mh) 2 



(52) 



The self-energy is due to the Hartree interaction between 
electrons of opposite spins; the Hartree and Fock inter- 
action between equal electrons cancel each other. 



The Dyson equation ( 50 ) involves the non-interacting 
Green's function Eq. ( 20 ) transformed to fc-space which is 



diagonal in the momentum and lead indices, G^j. p q (uj) = 
2ir6(k - q)6 a pg ( °l{w) with 0«>(w] given by Eq. @ 



We solve the Dyson equation (50) for G aka , k ,(uj) in 
fc-space and transform the result back to real space. 
The real-space Green's function G s ^ (x , x' , u>) then as- 
sumes the same form as for t he n on-interacting case (cf. 

(x, x', u>) as given by Eq. ( 20 )) but with renormalized 
scattering coefficients 



tak — t a (k-U/hVF)i 



r a k 



a(k-U/hv F )' 



(53) 
(54) 



This renormalization of the scattering amplitudes then 
corresponds to a mean-field sh ift of th e re sonance by 
U/Hvf', the general results Eqs. (27) and (31) correspond 
to the first order expansion in Uq of this result. Higher- 
order terms in the self-energy may not preserve the struc- 
ture of the non-interacting Green's function, in which 
case they cannot be cast into a simple renormalization of 
the scattering coefficients. 

The single-particle Green's function calculated above 
provides us with the corrections to the current j and 
the reducible current-current correlator Saa d \ cf. Eqs. 
( 36 1 and ( 39 1 . The non- interacting results for these two 



quantities, Eqs. (37) and (41 ), are modified by substitut 



ing the transmission probability T k with the expression 
Tk = T/ k _fj /nv F ) shifted in energy due to the interaction. 
For zero temperature (i? = 0) and in the linear-voltage 
regime, we obtain 



J |#=Q,V->-0 — -J- J-kp V) 



(55) 



S£ d)S< 'U=o,^o = j(l~f kF )f kF \eV\. (56) 

The irreducible current-current correlator Sa™ is de- 
termined by the two-particle Green's function. Here, we 



stay with the simple first-order correction Eq. ( 43 ) , as a 



proper resummation involves more complex diagrams in- 
cluding vertex corrections. The irreducible contribution 
to the current-current correlator is limited to opposite 
spins (where the reducible part does not contribute) and 
we obtain the result 



2 4C/ 
dk\dk,2 1 



s ti(°)l*=o = e h v<yl ~ 71 } sign(/iL ~ Mr) 

fen 

X 



(57) 



Afci/7 



(2tt) 2 7 2 (1 + (Afc!/ 7 ) 2 ) 2 (1 + (Afc 2 / 7 ) 2 ) 2 ' 



which for small voltage reduces to 
,4K 



S£(0)| 



i?=0.V->0 



°-v(l-V 2 ) 2 (f-) sign(F) 



(l + (Afc F ) 2 ) 4 



(58) 



Its sign 
sign(S£(0)| tf= 



■o,v- 



J = sign(77)sign(y)sign(Afc F ) (59) 



can be understood on the basis of the mechanism de- 
scribed at the end of Sec. IIIB2 The result (56) pro- 



vides the lowest order in V and Uq result for the equal- 
spin noise correlator; the opposite-spin correlator is gen- 
erated by the combination of interaction and asymmetry 
and is given by Eq. (57). 



It is instructive to compare our work with the results of 
Gogolin and KomnilP^. These authors calculated the log- 
arithm of the generating function log \ f° r the full count- 
ing statistics (FCS) in the Anderson impurity model to 
second order in the interaction parameter U /Y (with Y 
the width of the level), using a particular diagrammatic 
technique, in linear response (V — > 0) with the following 
results: i) The FCS assumes the same form as for non- 
interacting electrons but with renormalized transmission 
probabilities, ii) In the non-linear regime (finite voltage 
V) , spin-pairing correlations are generated to order V 2 in 
bias and to order (U/Y) 2 in interaction (pairing correla- 
tions have been founcP^ as well for a dot described within 
a master equation approach valid at large interaction U) . 
These results are obtained for a symmetric dot. In our 
work, we extend the analysis to an arbitrary scattcrcr 
with capacitive-like Coulomb interaction in the scatter- 
ing region. Rather than FCS, we limit ourselves to the 
calculation of the first two moments (current and noise) 
of the FCS and to first order effects in interaction. We 
confirm the renormalization of the single-particle scat- 
tering matrix in the linear response regime. For finite 
voltage V, we find that pairing correlations for electrons 
with opposite spins are enhanced for an asymmetric scat- 
terer, i.e., a scatterer with different reflection coefficients 
fLk 7^ r Rk for electrons incident from the left (L) and 
from the right (R) lead, as they appear already in the 
first order in interaction U/Y. 
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V. CONCLUSION 

We have extended the scattering matrix approach to 
include effects of interaction; while no restrictions have 
been imposed on the scatterer, our formalism applies to 
the special case of an interaction which is constant over 
the scattering region. Under these conditions, the fully 
interacting Hamiltonian can be expressed through the 
non-interacting (Lippmann-Schwinger) scattering states 
and scattering coefficients, hence the scattering aspect of 
the problem is dealt with in an exact way. On the part of 
the interaction, no knowledge is required about the wave 
function within the interaction region. The approach can 
be extended to the case where scattered electrons inter- 
act with electromagnetic and crystal degrees of freedom, 
photons, plasmons, and phonons. 

In our analysis above, we have considered a constant 
interaction kernel and one may ask, how well this choice 
approximates the situation when the kernel is short- 
ranged. For the single resonance level discussed in Sec. 
WJ the shape of the interaction kernel is not relevant; 



the matrix element Eq. ( 45 ) derives from the expression 



Eq. (151 for the finite scatterer in the limit of vanish- 



ing extent d and any interaction kernel is always long- 
ranged. For a finite dot of dimension d, the shape of 
the interaction kernel U (x, x') becomes relevant when its 
width drops below d. Here, we discuss the situation with 
sharp resonances and wave functions close to those of an 
isolated dot. As long as only a single resonance is in- 
volved, the modifications can be captured by a simple 
renormalization of the interaction strength Uq, with the 
Hartree and Fock contributions of similar weight. The 
Hartree contribution remains rather robust when includ- 
ing other occupied resonances, with a uniform renormal- 
ization of the interaction strength over the resonances 
(for a narrow kernel, the contributions from the low- 
est resonances, where the width of the kernel resolves 
the structure of the wave functions, is modified). For 
the constant kernel, the inter-resonance Fock terms are 
strongly suppressed as compared to the Hartree term, 
a consequence of the approximate orthogonality of the 
wave functions; a narrow kernel then enhances the rela- 
tive importance of these (actually small) inter-resonance 
Fock terms. Overall, changing the kernel, we expect 
small effects on the renormalization of the scattering co- 
efficients, on the transport current, and on the reducible 
part of the zero-frequency noise. For these quantities, 
the relevant physics involves electrons passing through 
one resonance and interacting with all the electrons on 
occupied resonances, which can be roughly included into 
a renormalization of Uq- Similarly, a small effect is ex- 
pected in the irreducible part of the zero-frequency noise 
as long as only one resonance is involved in the trans- 



port. However, at larger voltage, more than one reso- 
nance contributes to the irreducible noise correlator and 
the interaction between electrons passing through dif- 
ferent resonances becomes important; in this case, the 
specific enhancement of the inter-resonance Fock terms 
should be accounted for, e.g., via the introduction of ad- 
ditional 'cross-capacitances'. For an open dot with wide 
resonances and a short-range kernel, the wave functions 
in the dot have to be found explicitly when calculating 
the overlap matrix elements ( 12 ), the number of parame- 



ters increases, and the results can no longer be expressed 
through the scattering coefficients alone in the simple 
way we have done above. 

As a first application of our approach, we have deter- 
mined the one- and two-particle Green's functions within 
first order perturbation theory in the interaction param- 
eter Uq for a general scatterer - an analysis to higher or- 
der is more easily carried out for special scatterers, e.g., a 
single resonance level where the interaction Hamiltonian 
factorizes in fc-space. While the linear-response correc- 
tions to the one- and the reducible part of the two-particle 
Green's functions can be cast into a renormalization of 
the scattering coefficients, additional terms of different 
form show up in the irreducible part of the two-particle 
Green's function. In the non-linear regime at high volt- 
age, corrections appear which are beyond renormaliza- 
tion of the scattering matrix. Similarly, the interaction- 
induced corrections to the current and reducible part of 
the noise are accounted for through (the same) renormal- 
ized scattering coefficients. In addition, the interaction 
generates an irreducible contribution to the noise cor- 
relator which goes beyond renormalization of scattering 
corefficients. The latter gives rise to opposite-spin cur- 
rent correlations already to first order in Uq for the case 
of an asymmetric scatterer, an effect which has not been 
addressed so far. 

Our results are in agreement with those of Gogolin and 
Komnik 20 who reported a similar renormalization (in lin- 
ear response and to second order in Uq) of scattering 
coefficients in their study of the full counting statistics 
of the (symmetric) Anderson model; their spin-pair cor- 
relation appears only in second order of the interaction 
as their scatterer is symmetric. While our study is first 
order in Uq, it applies to any scatterer with arbitrary 
scattering coefficients. We hope, that the simplifications 
which come with our approach can be used to generate 
new (as with the asymmetric scatterer) or more precise 
(as with the resummation of a class of diagrams for the 
single resonance level) results in other situations. 

We thank Andreas Komnik for illuminating discus- 
sions and acknowledge financial support from the Swiss 
National Foundation, the Pauli Center for Theoretical 
Physics, and from the RFBR Grant # 11-02-00744-a. 
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